Hidden correlations of exciton complexes in self-assembled quantum dots and rings 
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The binding energies of trions (X + , X~) and biexciton (XX) in self-assembled semiconductor 
quantum dots (QDs) are very sensitive to the geometry and chemical composition of the QDs, and 
are random from dots to dots. However, in this letter, we show through analytical and numerical 
methods that the transition energies of the exciton complexes in self-assembled quantum dots and 
rings follow a simple and robust rule, i.e., the sum of exciton and biexciton transition energies 
minus the transition energies of trions is always positive and almost a constant for the same type 
of quantum dots and rings as a consequence of a pure Coulomb correlation effect. More interest- 
ingly, this quantity show a sharp transition when the topology change from a dot to a ring. This 
hidden correlation effect, directly measurable in experiments, offers a useful way to understand the 
photoluminescence spectra of self-assembled quantum dots and rings. 

PACS numbers: 78.67.Hc, 73.21. La, 68.65.Hb 
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New physics emerge when the size of a solid system 
reduces to nano-scale. In the self-assembled semiconduc- 
tor quantum dots (QDs), the three-dimensional confine- 
ment effects lead to atom-like electronic structure and 
long- living coherent quantum states. The confinement 
also enhances the Coulomb interactions among the quasi- 
particles, leading to novel physics in QDs other than bulk 
materials, e.g., Coulomb blockade effects (lj, and the non- 
Aufubau filling order for holes [2|, |3j in the InAs/GaAs 
quantum dots. The unique properties of QDs are not only 
of special interests in the view of fundamental physics, 
but also have important applications in quantum infor- 
mation processes. In these applications, the exciton com- 
plexes, including exciton (A), biexciton (AA) and trions 
(A + , X~) play extremely important roles. For example, 
neutral excitons can be used to generate single-photons 
[1, [f|, whereas biexcitons can be used to generate en- 
tangled photon pairs @, @] ■ Trions can be used to write 
in/read out the information of spin qubit, or to manip- 
ulate the spin states [1, 0, [l^, [ll[ . Due to the enhanced 
Coulomb interactions, the transition energies of the biex- 
citons and trions have significant (a few meV) energy 
shifts (i.e., binding energy) relative to those of the neu- 
tral excitons. Question arise that if we can find some 
simple relations between the transition energies of these 
exciton complexes? 

Unfortunately, both experimental measurements 12|, 
131 ] and theoretical calculations 14| show that the bind- 
ing energies of the exciton complexes change dramati- 
cally, even their signs, with respect to the sizes, shapes 
and chemical compositions of the QDs. Even worse, the 
binding energies seems to be random from dot to dot. 
It seems hopeless to find some simple relations between 
the transition energies of the exciton complexes. In this 



letter, we show however, through analytical analysis and 
numerical calculations that there is indeed a simple and 
robust relations among the transition energies of the exci- 
ton complexes in self- assembled quantum dots and rings. 

The binding energies of exciton complexes vary mainly 
due to the competition of the direct electron-electron, 
hole-hole and electron-hole Coulomb interactions (lil. fH] . 
which strongly depend on the size or the composition 
of a QD. To eliminates the influence of direct Coulomb 
interactions, we define a quantity A, 



A = A + AA — X + - X 



(1) 
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where, A, AA, A + , and A~ are the transition ener- 
gies of excitons, biexcitons, positive and negative trions, 
respectively. It is easy to show that A=0 under the 
Hartree-Fock (or single configuration) approximation, as 
the direct Coulomb energies and exchange energies cancel 
each other for the four types of exciton complexes. How- 
ever, A is generally not zero when we include correlation 
energies. 

We have calculated the photoluminescence (PL) spec- 
tra for a large amount of QDs with different sizes 
(radii and heights), compositions (InGaAs/GaAs or In- 
PAs/InP) and shapes (lens, cone, elongated and pyrami- 
dal dots). We first obtain the single particle energy lev- 
els and wave functions of the geometry-optimized QDs 
using an empirical pseudopotential method (IBj . where 
the total pseudopotential of the system are superpo- 
sition of the local, screened atomic pseudopotential of 
all (dot+matrix) atoms and the nonlocal spin-orbit po- 
tentials. The pseudopotential Schrodingcr equation is 
solved via the linear combination of bulk bands (LCBB) 
method [l6| . Many-body effects are included via the con- 
figuration interaction (CI) method [l7| by expanding the 
total wavefunction in Slater determinants for single and 
bi-excitons formed from all of the confined single-particle 
electron and hole states. 
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FIG. 1: (color online) (a) Calculated A for InGaAs/GaAs 
QDs with various sizes and geometries, (b) Experimental A 
extracted from literatures, (a), (▼) ( •), (♦), and (■) are 
extracted from Ref . [l2l Il3l . Ha . figL and [20l . respectively. 

Figure QJa) depicts the calculated A in InGaAs/GaAs 
QDs with respect to their exciton energies. It clearly 
shows that A, in contrast to the binding energies, are 
always positive and distributed in a narrow energy range 
from 0.8 to 1.8 meV even when the exciton energy change 
dramatically from 0.9 to 1.4 eV, with an average value 
equal to about 1.22 meV. We see that A hardly depends 
on the size, composition, and even the shape of the QDs. 
Therefore A represent a simple and robust relationship 
between the transition energies of exciton complexes. To 
verify the calculations, we show in Fig. [ljb) experimen- 
tal values of A extracted from the InGaAs/GaAs QDs 
PL spectra in the literatures [H O, 0, [H, H3| of four 
different groups. Indeed, the experimental A distributes 
in 1.0 - 2.0 meV as the exciton energy varies from 0.9 
to 1.4 eV, which strongly support the theoretical predic- 
tions. We have also calculated the A for InAs/InP QDs 
and find that A are distributed in 0.6 - 1.3 meV. 

Vanishes in the HF approximation, A is therefore a 
pure correlation effect. It not only represents a simple 
relation between the transition energies of the excition 
complexes, but also provides an easy way to character- 
ize the correlation effects in different nano-structures, 
which was previously only available in theoretical cal- 
culations. In order to gain an analytic understanding of 
these results, we take use of a perturbation theory. We 
use a typical lens-shaped InAs/GaAs QD as an exam- 
ple. Other kind of dots can be studied in a similar way. 
Since A=0 under Hartree-Fock approximation, A can be 
alternatively calculated as, 

A = A xx - A x + - A x - , (2) 



FIG. 2: A schematic figure of the excited configurations that 
couple substantially to the ground state configurations for 
XX, X + and X~ , divided in three categories (see text for 
detail). The numbers of percentage in each configuration is 
the weight of this configuration in the total many-particle 
wavefunctions from CI calculations. We omit irrelevant single 
particle levels. 

where, , Ax- , and Axx are the correlation energies 
(i.e., the energy difference between CI and HF calcula- 
tions) for the trions and biexcitions, respectively. For 
example, A x + = E C i(X + ) - Euf(X + ). To get the cor- 
relation energies for each type of excitons, we first solve 
the many-particle Hamiltonian in a single configuration 
approximation and obtain the energy Ej for j-th con- 
figuration To the second order approximation, the 
correlation energy of an exciton is, 

x ^ |(<j> (AQ|g/|^PO)l 2 ™ 

where, Hi is the Coulomb interactions among the quasi- 
particles. \<f>o(X)) and £^ HF) (X) are the ground state 
configuration and its energy for an exciton. The correla- 
tion energies of W,and trions X + , X~ can be calculated 
in the exact same procedures. 

We show schematically in Fig. [5] the excited configu- 
rations that couple substantially to the ground states of 
XX, X + and X~. These configurations can be classi- 
fied into three categories. We define the configurations 
in which only one electron (or hole) is excited as type 
I configurations. Type II configurations are those that 
two electrons (or two holes) are simultaneously excited, 
whereas type III configurations contain one excited elec- 
tron and one excited hole. Configurations of more than 
two particle excitations do not couple to the ground state. 

Type I configurations can be further divided as elec- 
tron excitations (type I-e) and hole excitations (type I- 
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h). As we see from Fig. [2] for X + only type I-e exci- 
tations contribute significantly to the correlation energy, 
whereas for X~, only type I-h excitations contribute sig- 
nificantly. On the other hand, both I-e and I-h con- 
tribute to the correlation energies of biexcitons. The 
off-diagonal (OD) matrix elements of I-h excitation for 
XX is (Hi-h)xx = 2(e s ,h s ;e s ,h 2s ) + (h s ,h s ;h s ,h 2s ) 
and for X~ is (Hi- h ) x - = 2{e s ,h s ;e s ,h 2s ). Here, we 
used a short notation, (e s , h s ;e s , h 2s ) = 



^(»'i)0fc.(ra)0e.(r l )0fc a ,(r2)^ , 3 ... 

—d rid r 2 . (4) 



r 2 ) n - r 2 



where, <f> es (rx), <f>h,( r 2), 4>h 2s ( r 2) are the electron (e) and 
hole (h) single particle wavefunctions, and e(ri — r 2 ) is 
the dielectric function (l7j . For simplicity, we dropped 
the spin notation in the above equation, but included it in 
real calculations. Since electron and hole wavefunctions 
have similar dimensions in QDs jl5j |, (e s , h s ; e s , h 2s ) ~ - 
(h s , h s ; h s , h 2s ). Therefore (h s , h s ; h s , h 2s ) almost cancels 
(e s ,h s ;e s ,h 2s ). As a result, (Hj^h)xx is much smaller 
than (Hi-h)x- ■ This cancellation effect is also the rea- 
son that I-e excitations contribute little to the correlation 
energy of X + . For the other I-h configurations, the OD 
term vanishes in cylindrical dots , because (e s , /i s ; e a , h p ), 
(e s , h s \e s , hd)=Q due to the symmetry of the wavefunc- 
tions. 

The difference of the denominators of in calculating 
the correlation energies of XX and X~ [see Eq. (J3))] also 
contributes to A. The single-particle part of the denom- 
inators are the same for XX and X~ , and the difference 
comes from many-particle interactions. It turns out that 
for XX and X~, the difference between the two denom- 



inators is J; 



Mi 



J. 



hh 



-2s: 



i.e., the difference of Coulomb 



energies between the hole s state with itself and s with 
2s states. This energy difference is much smaller than 
the denominators themselves. The correlation energy of 
XX from I-h excitations has a additional factor of 2 due 
to spin degeneracy. Combining all above factors, we find 
that I-h contribution to the correlation energies |Aj^l | 
is much larger than |A^~^|. Using the numerically cal- 
culated Coulomb integrals and single particle levels, for 
a lens-shaped dot with base D=25 nm, height h=2 nm, 
A^~_ ~ 2.5 A x £. The I-h configurations' contribution 
to A is about 0.93 meV. 

We can apply similar analysis to type I-e configura- 
tions, and find they also make positive contribution to 
A. However, the contribution from I-e excitations are 
much smaller than those of I-h ones, because electrons 
have much larger level spacing than that of holes in 
the InAs/GaAs dots and therefore larger denominators 
(about 3 times) in Eq. [3] than those I-h excitations. For 
the same dot above, I-e configurations contribute 0.31 
meV to A. 

Type II configurations can also be sub-classified into 
Il-h and II-e excitations. The Il-h OD matrix elements 



for XX are (H n _ h )xx = (h S: h s ;hi,hj), where i, j 
are the indices of the excited single-particle levels. If 
i j, the configurations would not couple to the ground 
state configuration in cylidrical dots. For i=j configu- 
rations, the most contributions are from p orbitals ex- 
citions. The contributions from higher energy excitions 
are much smaller due to larger energy spacing. Interest- 
ingly, we find that X + has the same OD matrix elements 
for Il-h excitations i.e. (H° I D h ) x+ = (H° D h ) X x- As 
a result, only the difference between the denominators 
contributes to A. For type Il-h configurations, the differ- 
ence between two denominators ~ 2(J^ — Jjif). We find 
that AEff_ h {XX) is always larger than AEff_ h (X+), 
because the coulomb integrals always satisfy > Jj£. 
As a result, type Il-h configurations always make positive 
contribution to A, but usually much smaller than those 
of type I configurations. Similar analysis applies for type 
II-e configurations. The denominator AEff_ e (XX) > 
AEff_ e (X~). Therefore, type II-e configurations also 
make positive contribution to A. However, due to much 
larger single particle level spacing of electrons, the contri- 
bution from type II-e configurations is even smaller than 
that from type Il-h ones. 

Type III configurations contain one excited electron 
and one excited hole, as shown in Fig. [2] We find that 
the OD Hamiltonian matrix elements for for XX, X + 
and X~ are the same for each type of excitation, i.e., 

{H°°)xx = (H?g) x+ = (H?g) x - = (e.,A.;e i ,ft j >, 
where j, j are the indeces of the excited single particle 
levels in the configurations. The denominators in cal- 
culating the correlation energies for these three types of 
exciton complexes [see Eq. ([3|)] are also very close for 
type III configurations. After consider the spin degener- 
acy, the correlation energy of XX almost cancels those 
of trions, and type III configurations have very little con- 
tribution to A. 

For typical cylindrical dots, the hidden correlation A 
mostly comes from type I configurations. The type II 
contribution is less but significant, whereas the type III 
contribution is negligible. For example, for the lens- 
shaped dot with D—25 nm and h=2 nm, the perturbation 
theory gives A=1.51 meV, in which type I configurations 
contribute 1.23 meV and type II ones contribute the re- 
mainder 0.28 meV. As to the individual configurations, 
the type I hole 2s excitation have the most significant 
contributions to A. This non-trivial result suggests that 
not including the 2s orbital in the CI basis will cause huge 
error in calculating A. The perturbation theory gives a 
little smaller A than that from CI method (1.80 meV), 
due to ignoring higher order terms. We have applied the 
same analysis for other type dots, with different sizes and 
shapes (e.g. elongated dots and pyramidal dots). Even 
though the details of the contributions from each config- 
uration may be quite different, the final results of A are 
very close, which suggests that there are more fundamen- 
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FIG. 3: The correlation energy of the exciton Ax and the 
hidden correlation A as functions of Rin/Rout in a InAs/GaAs 
nanostructure. When Ri n / Rout <0, the nanostructure is a 
QD, and becomes a QR when Ri n /R Q ut >0. 



tal physics remain to be resolved in future studies. 

We now discuss another type of important semiconduc- 
tor nano-structures, quantum rings (QRs), which have 
non-simply connected topology. In a QR, the single- 
particle energy level is more sensitive to the inner radii 
of the ring, whereas the Coulomb interactions dependent 
more strongly on the outer radii of the rings [HI, l22j ]. 
The single-particle energy levels in a QR are much more 



condensed than in a QD [21|, |22j of similar sizes, and 



therefore the particles are expected to be more corre- 
lated. Figure [3] depicts the correlation energy of an exci- 
ton Ax and A with respect to the ratio of inner radius to 
outer radius (Ri n /Rout) in a InAs/GaAs QR. We adapt 
the ring model from Ref. 0. The outer radius of the 
ring is fixed at R ou t= 15 nm, whereas the inner radius 
Rin increase from -15 nm to 7 nm. The structure with 
Rin < is a QD and the structure with Ri n > is a QR. 
Figure [3] shows that when the structure changes from a 
QD to a QR, Ax increases dramatically from about 2 
meV to about 10 meV. Remarkably, almost a constant 
in the QD region, A shows a sharp transition when the 
topology changes from a dot to a ring, where A jumps 
from 1.8 meV in the QD region to 3.5 meV in the QR 
region and saturates when Rin/Rout exceeds 0.2. This ef- 
fect can be directly measured in a dot to ring transition 
experiments (23[. 

To conclude, we have introduced a hidden correlation 
function as the sum of exciton and biexciton transition 
energies minus the transition energies of trions for the 
exciton complexes in semiconductor nanostructures such 
as quantum dots and rings. Measurable in experiments, 
the hidden correlation provides a deep insight to the 
Coulomb correlation effects of the exciton complexes. We 
show that the hidden correlation energy is positive and 
almost a constant for the same type of quantum dots. Re- 
markably, it shows a sharp phase-transition-like behave 



when the topology changes from a dot to a ring. The 
fundamental physics that govern this intriguing behaver 
remain to be fully explored. 
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